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r*| ■ Abstract 

l — 'j This paper presents, for ternary lysozyme-Na2S04- water system, the 

' thermodynamic data extracted from the measured values of four ter- 

J> . nary diffusion coefficients and the Onsager reciprocal relations. The 

calculation for derivatives of solute chemical potentials with respect 
' to solute molar concentrations was made using the method presented 

C*~) ■ in pQ. This method is applicable to systems in which the molar con- 

centration of one solute is very small compared to that of the other, 
| like in our case. The approach is illustrated for the lysozyme chloride- 

"^3 ■ Na 2 S0 4 -water system at 25° C, pH 4.5 and at 0.6 mM (8.6 mg/mL) 

lysozyme chloride and 0.1, 0.25, 0.5, 0.65, and 0.8 M Na 2 S0 4 con- 
centrations. The calculated solute chemical potential derivatives were 
pH * 1 used to compute the protein cation charge approximately. We also 

Q_|. compute the diffusion Onsager coefficients (Lij) for each composition 

j>! \ at pH 4.5. 
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Motivation 



In order to determine the protein structure through X-ray diffraction, high 
quality crystals are required. The protein crystallization process usually oc- 
curs in aqueous solution that contain salts as precipitant agents. The phys- 
ical and chemical properties of these solutions affect drastically the nucle- 
ation and crystal growth processes. Protein aggregation depends on protein- 
protein and protein-precipitant interactions in the solution. In these inter- 
actions, the effective charge of the protein plays an important role. From 
titration experiments [2] only the stoichiometric value of the protein charge 
can be determine, but it does not take into account the presence of ions 
that may bind on the macromolecules and change their net charge as already 
shown by diffusion experiments i3j. 

The diffusion of protein is one of the fundamental processes occurring in 
biological systems, and it is also an important step in the crystallization 
mechanism. Crystallization is an intrinsically non-equilibrium process, and 
concentration gradients will occur around the crystal. The protein crystal- 
lizes, reducing its concentration at the moving face of the growing crystal 
and creating a protein gradient between the bulk solution and the crystal. 
This gradient in turn causes multicomponent diffusive transport of protein 
and precipitant. Diffusion in protein crystal growth inevitably occurs under 
conditions for which no species has an uniform concentration raising the issue 
of multicomponent diffusion. 

The complete description of an n-solute system requires an n x n matrix 
of diffusion coefficients relating the flux of each solute component to the 
gradients of all solute components Q|. The importance of other species on 
protein diffusion follows from the one-dimensional flux relations [3]: 



in which the cross-term diffusion coefficients (off-diagonal elements (Dij) v i 7^ 
j) can be positive or negative. In ternary systems (n = 2), our case, the one- 
dimensional flux relations could be written as: 
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where J\ and J 2 - the protein flux and respectively salt flux, (Du) v and (D 2 2) v 
- the main-term diffusion coefficients relating to the flux of component to its 
own concentration gradient, and (D 12 ) v and (D 2 i) v - the cross-term diffusion 
coefficients relating the flux of each component to the gradient of the other. 
The index v from the diffusion coefficients shows that the experiment were 
done under the assumption the volume change on mixing and changes in 
concentrations across the diffusion boundary were small. Consequently, with 
a good approximation, the measured diffusion coefficients may be considered 
to be for the volume-fixed reference frame [3J defined by: 



where Vi is the partial molar volume of the ith species, and the subscript 
denotes the solvent. 

The importance of multicomponent diffusion has been recognized in the crys- 
tal growth community [SJ U\ and a crystal growth model has properly ac- 
counted for multicomponent diffusive transport in lysozyme chloride-NaCl- 
water system [SJ . The experimental multicomponent diffusion coefficients 
are essential for accurate modeling of protein transport, especially in view 
of the very large cross-term coefficient (-D 21 )„ reported here. Moreover, the 
concentration of supporting electrolyte dependence of all the diffusion coeffi- 
cients should be important for supersaturation region and also for its directly 
contribution to the protein flux. 

The use of ternary diffusion coefficients to determine binding coefficients and 
other thermodynamic data is very well established [HI EI] • For our ternary 
system, the molar concentration of one solute is very small compared to that 
of the other, and also small enough that an inverse concentration dependence 
dominates certain activity coefficient derivatives. For such systems, using 
the Onsager reciprocal relations (ORR), along with precision measurements 
of ternary diffusion coefficients from our earlier paper we determined 
concentration derivatives of the chemical potentials = dfii/dCj) of two 
solutes with respect to one solute molar concentrations. 
In our ternary experiments, the molarity of lysozyme chloride is small com- 
pared to that of Na 2 S0 4 . Thus, the self-derivative for lysozyme chloride, /in, 
is dominated by its concentration term. The self-derivative /X22 for Na2S04 
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is essentially that of the binary with minor correction. In order to obtain the 
molarity cross-derivatives (/112 and \i 2 \ which are unequal ^2]) we have to 
use two additional equations: (1) equality of the molality cross-derivatives 

dm 2 dmi 

where mj is the molality of solute i, and (2) the ORR equation: 

l^ii(D 12 )o - fiiz{D 2 i)o = ^22{D 2 i) - fi 2 i(D 22 ) (6) 

relating the four molarity derivatives and the ternary difussion coefficients in 
a solvent-fixed reference frame (Dij) Q [T4*l I15j . This method [I] yields to an 
estimate of lysozyme charge and a functional approximation to the change 
of the chemical potential of lysozyme chloride with Na2SC>4 concentration. 
This, together with the diffusion coefficients, will permit the modeling of 
protein crystallization processes. 

Experimental section 

All the experimental work was performed al Texas Christian University, in 
the Chemistry Department. 

Materials. All the materials, solution preparation procedures, apparatus 
and density measurement procedures are described in the work jH]. We used 
a hen egg-white lysozyme, recrystallized six times purchased from Seikagaku 
America. 

The molecular mass of the lysozyme solute, M\, was taken as 14307 g/mol, 
and this value (THj was used to calculate all concentrations after correction 
for the moisture and chloride content. Buoyancy corrections were made with 
the commonly used lysozyme crystal density [T71ITB1HHI of 1.305 g/cm 3 . 
The molecular mass of water, M , was taken as 18.015 g/cm 3 and the molec- 
ular mass of Na2S04, M 2 , was taken as 142.037 g/mol. 
Mallinckrodt reagent HC1 (~ 12 M) was diluted by half with pure water 
and distilled at the constant boiling composition. This resulting HC1 solu- 
tion (~ 6 M) was then diluted (pH 1.2) and used to adjust the pH of solution. 
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Preparation of Solutions. All solutions were prepared by mass with ap- 
propriate buoyancy corrections. All weighings were performed with a Met- 
tler Toledo AT400 electrobalance. Since the as-received lysozyme powder 
was very hygroscopic, all manipulations in which water absorption might be 
critical were performed in a dry glove box. Stock solutions of lysozyme were 
made by adding as-received protein to a pre-weighted bottle that had con- 
tained dry box air, capping the bottle, and reweighing to get the weight and 
thus mass of lysozyme. Water was added to dissolve the lysozyme, and the 
solution was weighed. An accurate density measurement was made and used 
to obtain the molarity of the stock solution. 

For ternary experiments, precise masses of Na2SC>4 were added to flasks con- 
taining previously weighed quantities of lysozyme stock solutions. These 
solutions were mixed and diluted to within 10 cm 3 of the final volume. The 
pH was adjusted, and the solutions were diluted to their final mass. 

Measurements of pH. The pH measurements were made using a Corn- 
ing model 130 pH meter with an Orion model 8102 combination ROSS pH 
electrode. The meter was calibrated with standard pH 7 and pH 4 buffers 
and checked against a pH 5 standard buffer. 

Density Measurements. All density measurements were made with a 
Mettler-Paar DMA40 density meter, with an RS-232 output to a Apple n+. 
By time averaging the output, a precision of 0.00001 g/cm 3 or better could 
be achieved. The temperature of the vibrating tube in the density meter was 
controlled with water from a large well-regulated water bath whose temper- 
ature was 25.00± 0.01 ° C. 

Free-Diffusion Measurements. For binary Na2S04-water and ternary 
Lys-Na2S04- Water we performed measurements for free-diffusion using the 
high-precision Gosting diffusiometer HH 122] operated in its Rayleigh in- 
terferometric optical mode. The procedure for measuring binary (D2)v and 
ternary diffusion coefficients (Djj)„ were described in detail in the work [8]. 
In order to measure the four diffusion coefficients of the system, the exper- 
iments must be performed with at least two different concentration differ- 
ences at each combination of mean concentration [201 123 121] ■ For ternary 
experiments, for each pair of mean concentrations, two measurements were 
performed with a± — and the two with a\ = 0.8 (a^ - the refractive index 
fraction due to the ith solute [8J). 
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In order to make the data analysis of the free-diffusion experiments we used 
the Fick's second law: 



dt 



2 



d 2 d 
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dx 2 



i = 
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(7) 



for two solutes. We made the assumption that the concentration differences 
of the solutes across the initial boundary are small enough and the diffusion 
coefficients are constant Also the volume changes on mixing were neg- 
ligible, thus all the measured diffusion coefficients are given relative to the 
volume-fixed frame of reference defined by equation (4). Also we could made 
a comparison with the dynamic light scattering [27| results from our previous 
paper [IT]. 



Ternary diffusion experiments were performed on the lysozyme chloride- 
Na2SC>4-water system at pH= 4.5 and T=25° C. The four ternary diffusion 
coefficients obtained, (Du)„, (D 2 2)i>, (Di2)« an d (D 2 i)„ were published in our 
earlier paper [TT] for a mean Na 2 S0 4 concentration of 0.1, 0.25, 0.5, 0.65 and 
0.8 M. The interferometric data for the diffusion coefficients (D n )„, (D22)j;, 
(Di 2 )„ and (D 2 i)„ are reported in the Tables 1,2,3. 

Partial molar volumes values, V\, V 2 and V a , were calculated for each 
component using eqs A- 7 (q = 2) and 5 in [T3j and reported in the Tables 



Values of mean density d and Hi = (dd/dCi)T, p ,c4,j^i m the Tables 1,2,3 
were calculated using densities of all eight solutions from each experiment 
set. Densities were assumed to be linear in solute concentrations respecting 
the equation [S]: 



where C\ and C2 are the averages of the mean concentrations for all four 
experiments in a series. 

Use of irreversible thermodynamics and diffusion data to calcu- 
late chemical potential derivatives. For our ternary system lysozyme 



Results 



1,2,3. 



d=d+H 1 (C 1 - C x ) + H 2 (C 2 - C a ) 



(8) 
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chloride-Na2S04, the molar concentration of a macromolecule and the sup- 
porting electrolyte are small and large, respectively. So, in this case it is 
possible to estimate the chemical potential derivatives fin and [i 22 [I]- 
The molality cross-derivative relation, eq.(5), comes from classical thermo- 
dynamics. From eq. (5), an expression can be derived relating the four mo- 
larity partial derivatives /ijj The eq.(6), the ternary ORR of irreversible 
thermodynamics, relates the four (Dij) values of diffusion coefficients in a 
solvent-fixed reference frame and the four fiij values . 



Fundamental equations. The analysis of chemical potential derivatives 
was made in terms of quantities referred to a solvent-fixed reference frame, 
identify by a subscript 0, with the diffusion Onsager coefficients denoted 
by (Lij) . The solvent-fixed (Dij) a values shown in the Tables 1,2,3 were 
obtained from experimental volume- fixed {Dij) v values by standard equations 
involving the % [U EH EH] . 

The diffusion Onsager coefficients (Lij) and the solvent-fixed (£)y) diffusion 
coefficients are related (in matrix form) by: 



(D u )o (D 12 ) 
(D 21 ) (D 22 ) 

(£11)0/^11 + (-^12)oAt21 {L\l)o^l2 + (Li 2 ) fi 22 

(L21 )oA f n + (L 22 ) fi 2 i (L 2 i) fi 2 i + (L 22 ) Q fi 22 



(9) 



The inverse realtion is: 



1 



(-^11)0 (£12)0 

(£21)0 (£22)0 J /ill/i 2 2 - ^12/^21 

yU22(-Dn)o - A i 2i(-Di2)o Vn(D 12 ) - /ii 2 (-Dn)o 
^22(D 21 ) - i~l 2 i(D 22 ) jj n (D 22 ) - fi 12 (D 21 ) 



X 



(10) 



Since the ORR, (L 12 ) Q = (L 2 1) Q [T3J[T3] apply to the solvent-fixed frame, eq. 
(10) yields eq. (6). 

Also, from eq. (5) for cross-derivative molality and the relations between Cj 
and rrii, we can show that [TTj : 



dfii 



C M 8m 2 

ii2i0--c 1 v 1 



/ii 2 (l - c 2 v 2 ) 
fi 22 C 2 Vi 



\i xx G x V 2 = 

1 dfi 2 
C M dm l 



(11) 
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where M is the molecular mass of water. 

In order to calculate the /i^- molarity cross- derivatives, we took in account 
their general thermodynamics expressions in terms of volume concentrations 
and the corresponding mean ionic activity coefficients yi for volume concen- 
trations and also we assumed that lysozyme chloride has stoichiometry 
LyCL- p . Thus, the fiij for our case could be written in a matrix form as: 



/ill A*12 
A*21 A*22 



RT x 



'121 



j_ 

Ci 



z p /2C 1 +C 2 + ( Z P+ l 

z p |_ o 9 ln 



, dlnyx 

aci 



z p /2Ci+C 2 



dCi 



z p /2C 1 +C 2 
C 2 ^ z p /2Ci+C 2 



8C 2 
i Q d ln y 2 

^ z ac 2 



where the quantity z p /2C\ + C 2 is equivalent to the total normality N of our 
ternary solution. 

Using the eq.(12) we computed the partial derivatives of chemical potentials, 
/ijj, for the case in which the molarity of at least one component is very low. 
From eq.(12) we calculated the fin and /z 22 taking in account the terms which 
are dominant for our case. In order to calculate the other two derivatives 
chemical potentials, fi\ 2 and // 2 i we use the eqs. (6) and (11) and the four 
(Ai) c : 



A*12 



A*21 



{/in[CiF 2 (£> 22 ) - 
^22^2^1(^22)0- 
{(1-C 2 V 2 (D 22 ) 

{fin[CiV 2 (D 11 ) - 

// 22 [WiPii) - 
{(1-C 2 V 2 (D 22 ) 



l-dVx^D^),,]- 
(l-C.V^DM}/ 

-(l-CiVOpn)*;} 

l-C 2 V 2 )(D 12 ) ] - 

(i-c 2 y 2 )p 21 ) ]}/ 

-(1-CxVi) pn)o} 



(13) 



(14) 



All the values for derivative chemical potentials fJ>%j £1X6 reported in the Tables 
1,2,3 for each salt concentration. 

Using the eqs. (10), (12), (13) and (14) we calculated the thermodynamics 
transport coefficients (Lij) for pH 4.5 and C% = 0.60 mM, which are reported 
also in the Tables 1,2,3. 



Calculation of lysozyme chloride charge using the values for /ii 2 and 
yU 2 i is described in [I]. Looking into eq. (12), the expressions for /i 12 and \x 2 \ 
suggest that we can obtain the lysozyme chloride charge multiplying /i 12 /i?T 
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and by z p /2Ci + C 2 and we'll receive the following dependence: 



(z p /2Ci + C 2 ) 



^12 




(9 In |/i 



(15) 



dC 2 



Y 21 = (^/2d + C 2 ) 



At-21 



(16) 



The dependence of Y21 and Y21 on (z p /2C\ + C 2 ) is shown in the Fig.l. From 
the graph we could estimate the protein charge, z p and to compare it with 
the value obtain from thermodynamics data from pQ. At pH 4.5 ,the average 
value of z p obtained from eqs. (15) and (16) is z v = 3.29, in comparison 
with the value obtained, using the same approach, for lysozyme-NaCl- Water 
system, z p = 8.9. This data were obtained from thermodynamics data for 
ternary system, which were in turn obtained in part from transport data. 

Conclusions We reported the complete set of multicomponent diffusion co- 
efficients for ternary lys-Na2S04-water system at concentrations high enough 
to be relevant to crystallization studies, in the volume- fixed frame, (Dij) v and 
in the solvent-fixed frame, (Z?y) . Also we calculated the derivatives chem- 
ical potentials % for our ternary system and after that we estimated the 
lysozyme chloride charge, z p , from irreversible thermodynamics and diffu- 
sion data. We also reported the four thermodynamics transport coefficients 
(Lij) at pH 4.5 and G x = 0.60 mM. 
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Table 1 



Ci(mM) 


0.6000 


0.6000 


C 2 M) 

(/(gem -3 ) 

^(lO^mo/" 1 


0.1000 


0.2500 


1.012189 


1.030671 


4.286 


1.149 


H 2 (10 3 gmol" 1 


0.12500 


0.12230 


VAcirfimol -1 ) 


10050 


10182 


V 2 {cm 3 mol~ l ) 


17.090 


19.780 


V {cm 3 mol~ 1 ) 


18.067 


18.058 


(DnUlO-Ws- 1 ) 


0.1169± 0.0001 


0.1090± 0.0001 




-0.000013± 0.000001 


0.000108± 0.0000 


(D 21 ) v (10^m 2 s- 1 ) 


2.49± 0.01 


4.14± 0.01 


(D 22 ) v (lO- 9 m 2 s- 1 ) 


0.9661± 0.0001 


0.8826± 0.0001 


(DnUlO-^s- 1 ) 


0.1176 


0.1097 


(D 12 ) o (10- 9 m 2 s- 1 ) 


-0.0000031 


0.0001193 


(DMIO-Ws- 1 ) 


2.6 


4.5 


(L> 22 ) o (10- 9 m 2 S - 1 ) 


0.968 


0.887 


Vn/RTiM- 1 ) 


1787 


1715 


Vn/RTiM- 1 ) 


41.9 


22.5 


fi^/RTiM- 1 ) 


64.2 


43.3 


fi^/RTiM- 1 ) 


21.894 


8.131 


RT( y L 11 )o(10- 9 Mm 2 S - 1 ) 


0.0000707 


0.0000684 


RT(L 12 ) o (l0~ 9 Mm 2 s- 1 ) 


-0.0001354 


-0.0017452 


RT(L 22 ) o (10- 9 Mm 2 s- 1 ) 


0.044487 


0.1094762 
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Table 2 



Ci(mM) 


0.6000 


0.6000 


C 2 (M) 
<i(gcm~ 3 ) 


0.5000 


0.6500 


1.060538 


1.078032 


ifi(lO 3 0moZ _1 


4.104 


4.049 


fi 2 (10 3 omoZ- 1 


0.11733 


0.11610 


V 1 (cm 3 mor 1 ) 


10209 


10257 


V2(cm 3 mo/ _1 ) 


24.720 


25.930 


y o (cm 3 mol _1 ) 


18.026 


18.013 


(AO^io-Ws- 1 ) 


0.0969± 0.0001 


0.0894± 0.0001 


(L» 12 ),(10- 9 m 2 S - 1 ) 


0.000132± 0.000001 


0.000134± 0.000001 


(D 21 ) v (10- 9 m 2 s- 1 ) 


6.75± 0.01 


8.26 ± 0.01 


(D^UlO-Ws- 1 ) 
(DnUW-Ws- 1 ) 


0.7791 ± 0.0001 


0.7294 ± 0.0001 


0.0977 


0.0902 


(D 12 ) o (10- 9 m 2 s- 1 ) 


0.0001444 


0.0001475 


(D 21 ) o (10- 9 m 2 s- 1 ) 


7.3 


9.0 


(L> 22 ) o (10- 9 m 2 S - 1 ) 


0.790 


0.742 


/i 11 / J RT(M- 1 ) 


1691 


1685 


li l2 /RT{M- 1 ) 


17.7 


17 




36.8 


35.6 




3.755 


2.811 


RT(L 11 ) o (10- 9 Mm 2 s- 1 ) 


0.0000634 


0.0000601 


RT(L 12 ) o (10- 9 Mm 2 s- 1 ) 


-0.0002606 


-0.0003110 


RT{L 22 ) o {l0- 9 Mm 2 s- 1 ) 


0.2117597 


0.31738664 
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Table 3 



Ci(mM) 


0.6000 


C 2 (M) 
d(gcm~~ 3 ) 


0.8000 


1.095745 


H^ltfgmol- 1 


3.954 


H 2 (10 3 gmol- 1 
Vi(cm 3 mol~ l ) 


0.11502 


10341 


V 2 (cm?mol~ l ) 


26.980 


V (cm 3 mol~ 1 ) 


17.993 


(DnUlO-Ws- 1 ) 
(DMIO-Ws- 1 ) 
(D 21 ) v (10^m 2 s- 1 ) 


0.0822± 0.0001 


0.000130± 0.000001 


9.50± 0.01 


(D 22 UlO~ 9 m 2 s- 1 ) 


0.6900± 0.0001 


(D 11 ) o (10- 9 m 2 s- 1 ) 


0.0829 


( J D 12 ) o (10- 9 m 2 s- 1 ) 


0.0001418 


(D 21 ) o (10- 9 m 2 s- 1 ) 
(L> 22 ) o (10- 9 m 2 S - 1 ) 


10.4 


0.706 


fin/RTiM- 1 ) 


1682 


li 12 /RT{M- 1 ) 


16.1 


^21/ RT(M~ 1 ) 


34.5 


fi^/RTiM- 1 ) 


2.241 


RT(L 11 ) o (10- 9 Mm 2 S - 1 ) 


0.1808868 


RT(L 12 ) o (10- 9 Mm 2 s- 1 ) 
RT(L 22 ) o (10- 9 Mm 2 s- r ) 


-0.000341 


0.3173864 
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